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Abstract 

We study a family of quadratic stochastic differential equations in 
the plane, motivated by applications to turbulent transport of heavy 
particles. Using Lyapunov functions, we find a critical parameter value 
ai = 0:2 such that when a2 > Oi\ the system is ergodic and when 0:2 < 
ai solutions are not defined for all times. Hormander's hypoellipticity 
theorem and geometric control theory are also utilized. 
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1 Introduction 

Understanding the ergodic behavior of both deterministic and random dy- 
namical systems is of paramount importance in applications. This is be- 
cause many natural phenomena are either in steady state or close to it; thus 
the long-time dynamics often reflects most accurately the situation at hand. 
Proving or disproving ergodicity in general, however, remains a challenge. 
In this work wc discuss useful methods, as applied to a specific family of 
stochastic differential equations. 



As in (7|[l7lll2 

, our main methodology is the use of Lyapunov functions. 
These are widely utilized in the literature, e.g. |5 14 21 , to study the behavior 
of diffusion processes. It is well known that finding a Lyapunov function to 
verify global stability can be highly nontrivial. We aim at a systematic 
approach to this problem by introducing the notion of a Lyapunov covering. 
This method consists of identifying the regions which the process must exit 
prior to leaving the state space. 

Moreover, we hope that our use of geometric control theory [I| [Io|jl3|[22] 
will elucidate the study of stochastic differential equations with degenerate 
noise. One class of such systems arises by writing a second-order equation 
in first-order form [6 ,14,20 . In such equations, noise is effectively trans- 



ferred through the average dynamics so that the support of the diffusion for 



all positive times is the whole space. In view of the support theorem 22 
this is proven by explicitly solving the associated control problem. In the 
cases treated here, it is more convenient to use geometric methods to un- 
cover the accessibility sets of the control system. We find that, due to the 
even (quadratic) order of the coefficients of the average dynamics, noise only 
spreads to certain strict subsets of the state space. 

In what follows, let ai, 02 G M, ai, ^2 > and ki >0,K2> 0. We study 
the stochastic differential equation in 



d2. 



dXt 



{a2Yt - 



dt + V2ki dBl 
a2XtYt) dt + y^dB^, 



(1) 



where B^.^'' and -Bp"* are independent standard Brownian motions. Consider 
a two-dimensional spatially smooth turbulent flow. It was noted in [2] that 
in our system with oi = 02 = —1 and ai = 1, 0:2 = 2, Xt and Yt model 
the transverse and longitudinal components of the velocity difference of two 
heavy particles transported by the flow. Ergodicity of this equation was 
assumed in [2] to extract information on particle clustering as expressed 
through the top Lyapunov exponent. Subsequently, ergodicity was proven 
in [gI. Here, among other results, we provide an alternate proof of this fact. 



It was also noted in [2||6] that, in this special case, equation (1.1) can be 
written succinctly in the complex variable Zt = Xt + lYt as: 



dZt 



'-Zt - Z^) dt + V2^i dBl 



(1.2) 



Here we show that ergodicity of equation (1.1) holds whenever 0:2 > ai] 
thus stability in (1.2) above is not due to its holomorphic structure. We 
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shall see that this result is optimal in the sense that when ai > a^-, there 
are solutions which reach infinity in finite time with positive probability. It 
should be noted that the choice ki > 0, K2 > is optimal as well. For fixed 
Oi, 02 G M and ai,a2 > 0, if /ti > 0, ^2 = 0, with probability one equation 



(1.1) has solutions which explode to infinity in finite time. This is a simple 
consequence of Feller's test applied to solutions starting along the negative 
real axis, hence we omit further discussion. 

By working in this generality, our goal is to make progress in under- 
standing turbulent transport of heavy particles in spatially rough fiows as 
well. In a two-dimensional fiow with spatial Holder exponent h G (0, 1), the 
components, Xi and F^, of the velocity difference now obey: 

dXt = {-Xt- R;\hX^ -Y^^))dt + y^idBl^^ (1.3) 
dYt = {-Yt~{l + h)R-^XtYt)dt + ^2dBf^ 
dRt = (1 - h)Rt dt, 

where Rt is proportional to |S't|^~'' where St is the particle separation at time 
t > [3]. For simplicity, to gain insight into the dynamics (1.3) we assume Rt 



is a positive constant; in which case, we see that equation ( |1.3[ ) falls within 



the class (1.1) 



The organization of the paper is as follows. In Section [2| we fix notation, 
state the main results, and outline our methods of proof. In Section |3} 
we show that when 0:2 > ai, the process {Xt,Yt) is nonexplosive and that 
equation (1.1) has (at least one) invariant probability measure. In Section |4| 



we prove uniqueness of the invariant probability measure using a combination 
of Hormander's theorem and geometric control theory arguments. In Section 
5| we prove that when ai > 02, the process {Xt,Yt) with initial conditions 
in a certain region goes to infinity in finite time with positive probability. 



2 Main Results 

Let us first fix notation and terminology. Throughout the paper B denotes 
the Borel a-field of subsets on M^. We assume the standard Brownian mo- 
tions, bI^'^ and bI'^\ are defined on a common probability triple {Q,J^,P). 
Zt denotes the two-dimensional process {Xt,Yt) defined by equation ( 1.1[ ). 



Zt is defined on (fi, J-", P) until the random time at which it leaves M . More 
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formally, for n G N, if we define 



r„ := inf{|^t| > n}, 

and let r be the finite or infinite limit of as n — )■ oo, is defined on 
{Q,J^,P) for all times t A t, t > (where here and in what follows a Ab 
denotes the minimum of a and b). To emphasize that both and r depend 
on Zq = z G M^, the superscript notation and is sometimes used [?]. 
We opt instead to account for the dependence on z in the measure via Pz 
and Ez- 

Using the notation above, Zt/^^ is a Markov process with transition kernel 

Pit,z,A) ■.= Pz{Zt;,,eA}, 

for t > 0, z G M^, and A E B. The measures P{t, z, ■ ) give rise to a semigroup 
{Pt}t>o which acts on real- valued bounded S-measurable functions $ via 

Pt<^{z) = / P{t,z,dw)<l>{w) 

and on finite i3-measures /x via 

fiPtiA)= I ii{dw)P{t,w,A), 



for A €z B. The fi-measure fi is called an invariant measure if for all A E B 
and t > 

fiPtiA) = f^iA). 

A positive invariant measure fi can be normalized to have total mass one; 
in which case, the resulting measure z/ is a probability measure that satisfies 
uPtlA) = z/(y4) for all A E B, t > 0. We call u an invariant probability 
measure. It is well known that invariant probability measures represent the 



long-time behavior of the process Zt^r- For a nice discussion of this, see 20 
Let 

d d d"^ d"^ 
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For a suitable class of functions, e.g. C^-functions $ : — )■ M, for n G N we 
have Dynkin's formula: 



[<^{Z,^,^)]-<^{z) = E^ 



tATn 



m{Zs)ds 



We call L the generator of the Markov process Zt^r- 

We now state our main results. The first says, provided 0^2 > ai, t = 
00 almost surely regardless of the initial condition ^ G M^. Moreover, the 
Markov process Zt-, which is now defined for all finite times t > and for all 

= -2 G M^, has a unique invariant probability measure. The second result 
shows that, in sharp contrast to the above, when 02 < cti, there exists a non- 
empty set A of initial conditions for which the process Z^ with Zq = z & A 
explodes with positive probability in finite time, i.e., Pz{t < 00} > for 
z E A. The set A will be explicitly described in Section |5| 

Theorem 2.1. 7/0:2 > ai, then: 

(1) For all z G M^^ 

Pz {r < 00} = 0. 
Thus for all z G M^^ P, {Zt = Zt^r for all t > 0} = 1. 

(2) There exists a unique invariant probability measure v. 

Theorem 2.2. If a\ > a2, then there exists 7^ A C such that for all 

z e A, 

P, {t <oo}> 0. 



Theorem 2A_ will be established by proving Lemma 2^ and Lemma 2.4 
below. Lemma 2.3 implies that Zt is nonexplosive and has at least one 



invariant probability measure (see [7,17,20]). Lemma 2.4, which is only nec- 
essary when Ki = 0, shows that the supports of any two ex^remaij^ invariant 
probability measures contain a non-empty open subset in common. This, 
coupled with regularity of the transition measures of the process Zt, proves 
uniqueness. For more details and further information, see [H[7|[8| [l5}jl7| . 



^An invariant probability measure v is extremal if whenever v ^ {1 
vi, V2 invariant probability measures and A G (0, 1), then vi = V2 ~ v. 



X)v\ + \v2 for 
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Lemma 2.3. If a2 > «i, there exists a C°° -function $ : — )■ M such that 

(1) $(2;) 00 as \z\ — )• 00. 

(2) There exist positive constants C, D such that 

L^{z) < -C^{z) + D, (2.2) 

for all 2 G M^. 
Lemma 2.4. 



p2 . 



(1) Fort>0, z e 

P(t, z, dw) = pit, z, w) dw, 

where p{t, z, w) is jointly C°° in (t, z, w) G (0, 00) x x and dw is 
Lebesgue measure on M^. Moreover, any invariant probability measure 
V (such measures exist by Lemma \2. 5| ) satisfies 

v{dw) = py{w) dw, 

where e C°°(M2). 

(2) Let fi,fi be extremal invariant probability measures. Then there exists a 

non-empty open U such that 

supp fi n supp ft D U. 



To establish Theorem 2.2, we use the following lemma, the proof of which 
is a simple consequence of Dynkin's formula. For more information, consult 

III- 

Lemma 2.5. There exist a bounded C°°-function \E' : [0,oo), strictly 

positive on A 7^ C M^, and a constant E > such that: 

L'^iz) > E^{z), 

for all zEM?. 
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We note that in the case when Ki > and 02 > "we can prove that 
the rate of convergence to the invariant probabihty measure u is exponential. 
To see this, first note that we may choose a constant F > such that 
= $ + F > 1. The rate can then be quantified by the norm || ■ ||$^ 
defined on ;B-measures /x via: 



sup 

0</<>I'f 



f{y) dfi{y) 



where the supremum is taken over nonnegative i3-measurable real-valued 
functions / bounded above by [8 , 17 . 



Corollary 2.6. If 1^1 > 0, then there exist r,C > such that 

\\P{t,z, ■)-uUp<C<l>F{z)e-'-\ 



for all t > 0, z & M.'^ . Moreover, since > 1 we can replace the left-hand 
side of (2.3) by \\P(t,z, ■) — i^\\tv where \\ 



(2.3) 

hand 

TV is the total variation norm 



on B-measures. 



Proof. The second conclusion easily follows from the first. To see the first, 
since both ki and K2 are strictly positive, it follows that 

supp P(t, z, ■)= 



for all t > 0, and all z G M^. Moreover, Lemma 2.4 (1) is satisfied. This 
together with Lemma 2.3 proves the result by Theorem A. 2 of Isl. □ 



3 Existence of a Lyapunov Function 



In this section, we prove Lemma |2.3[ Unfortunately, it is difficult to write 
down a single formula which would define $ for all z. There are a number of 
reasons for this. First, when ki = K2 = 0, there is an unstable trajectory in 



(1.1) along the negative X-axis. Thus some noise, e.g. K2 > 0, must be uti- 
lized so that solutions do not explode. In terms of estimating L$, this means 
second order terms must be used along this trajectory. Second, it is clear that 
noise with arbitrarily small magnitude in the F-direction will direct Zt with 
Zq = (Xo,lo) = (3^)0), a; < off the unstable trajectory. Once the process 
leaves the axis, the same noise that initially helped could now steer the pro- 
cess back to the unstable trajectory. This is reflected in the estimates for L$ 
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in a sufficiently large region which includes the negative X-axis. And third, 
once the process exits this region we must assure stability in the remaining 
regions as well. Because the nature of the unperturbed dynamics changes 
significantly from the previous region, this calls naturally for a function $ of 
a different type. We now introduce the following definitions to elucidate our 
procedure. 

Definition 3.1. Let U C be an unbounded region such that dU is a con- 
tinuous, piecewise smooth curve. We call a function (j) : U — )• M a Lyapunov 
function (for the operator L) on U if: 

(I) (j) G C^{U). 

(II) (f){z) — )■ oo as |2;| — )■ cxo, z E U. 

(Ill) There exist C,D>0 such that L(f){z) < -C(j){z) + D for all z eU. 

Definition 3.2. Let 0i, 02, • • • , 0Af be Lyapunov functions on Ui, U2, 
Un respectively and suppose there exists R> such that 

N 

i=i 

where is the open ball of radius R centered at the origin in M^. We call 
{(01, Ui), (02, U2), . . . , {4>N, Un)} a Lyapunov covering. 

To exhibit $, we will first find a Lyapunov covering. The existence of 



such a covering is of course not sufficient to prove Lemma [2. 3[ As suggested 
by the following proposition, however, finding a Lyapunov covering is a step 
in the right direction. 

Proposition 3.1. Suppose (f) is a Lyapunov function on U . Let V d U he 

another unbounded region, satisfying the same assumptions as U and such 
that dUndV = 0. // TV = mft>o{Zt G V^}, then for z eV 

PArv<r} = l. 

Proof. Let tu = inft>o{^t G U'^} and T^^uif) = A r^r A t. By adding a 
sufficiently large constant to 0, we may assure the resulting function is a non- 
negative Lyapunov function on U. Let us use to denote this non-negative 
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Figure 3.1: Cartoon of Proposition 3.1 




function. It is clear that the much weaker bound L(f){z) < C(j){z) + D is 
satisfied on U. Letting ip{z,t) = e~'"*(0 + ioi z E V we obtain by 
Dynkin's formula 



< 0. 



We see that 

Tn,u{t))] > e-^VnP. {Tn < Tu A t} , 

where (pn = inf^G(7n9B„ — )■ oo as n — )• oo. Combining the previous two 
estimates yields 

Pz{rn<ruAt}<(p-'e'''^{z,0). 

Letting ri — i- oo we see that Pz{t < tu A t} = for all t > 0, whence 
Pz{t > Tij} = 1. By path continuity of Zt, Pz{tu > Ty} = 1. Hence 
Pz{t > Ty} = I as claimed. □ 



Remark 1. Note that the same behavior is true if Property (III) of Definition 



3.1 is replaced by a much weaker bound. Indeed suppose U is as in Definition 
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3.1 and M satisfies (I) and (II) of Definition 3.1 If satisfies the 

bound 

L(f){z) < C<p{z) + D, 



for all 2; G for some C, D > 0, then the conclusion of Proposition 3A 
remains valid. 

The existence of a Lyapunov covering {{(t)j,Uj}^^i guarantees that Zt 
cannot go to infinity directly through any unbounded region V C Uj with 
continuous boundary such that dV fl dUj = 0. As we shall see in our case, 
it is possible to extract subsets Vj C Uj with continuous piecewise smooth 
boundaries such that dUj fl dVj = for all j = 1, 2, . . . , and Vj)}jL^ 
is a Lyapunov covering. 

Definition 3.3. We call any {{4>j, Vj)}jLi with the above properties a strong 
Lyapunov covering subordinate to {{(j)j,Uj)}jLi. 

Hence, given the existence of a strong Lyapunov covering V^OjjLi, 
Zt cannot go to infinity staying inside one of the Vj. Moreover 

N 

R^ = [j V, U Br, 

J=l 

for some R > 0. Thus it seems likely that Zf cannot diverge in finite time. 
However, it is possible that Z^ oscillates between two or more regions on 
its way to infinity in finite time. To eliminate this possibility we glue our 
Lyapunov covering together so that we have a globally-defined Lyapunov 
function $. 

The construction of $ is done in three stages: 
Stage 0: Reduction of parameters. 
Stage 1: Existence of a strong Lyapunov covering. 
Stage 2: Existence of a globally-defined Lyapunov function $. 

3.1 Stage 



Equation (1.1) depends on 01,02 G M, ai,a2 > 0, and ki > 0, ^2 > 0. 



We thus may expect our globally-defined Lyapunov function to depend on 
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these parameters as well. To simplify matters, let a = (ai, a2), a = (ai, 02), 
K = ^2), and L'^ct be the generator (2.1) of Zt that solves equation (1.1) 
with those parameter values. Define 

Af{a, CK, k) = 1$ : $ is a Lyapunov function on for L'^^.} ■ 

We have the following: 

Proposition 3.2. Let 02 > cti > and suppose there exists K2 > such 
that for all a e and ki > 0, Af{a, a, k) ^ 0. Then Mih, ol,l)^^ for all 
b G M^ t G M>o X M>o. 

Proof. Let t = {li,L2) G R>o x ]R>o and b G = (&i,&2) be arbitrary. 
Define c = (^262"'^) 3, a = cb, and k = {c^ii, ^2)- By assumption, there exists 
\l/ G J\f{a,a, k). Define $(2) = \l/(cz). It is easy to see that $ satisfies 



Properties (I) and (II) of Definition 3.1 in all of M^. To see (III), we use the 
chain rule to obtain: 



< 



for some C,D > 0. 



-C^{cz) + D 



□ 



Using Proposition |3.2[ it is enough to show, given 0:2 > «i, that there 
exists K2 > sufficiently small such that A/'(a,Q;,K) 7^ for all a G 
and Ki > 0. This is equivalent to saying that for all 0:2 > ai, there exists 
K,2 = i^2{a.) > and such that G A/'(a, a, k) for all a G and all 
Ki > 0. To obtain such a function $^q,, we first find a Lyapunov covering. 



3.2 Stage 1 

We will now construct families of strong Lyapunov coverings that depend on 
CK, a, and k. By Proposition 3.2, we will fix 02 > ai and pick K2 = ^2(0;) > 
such that {(pi, Ui) for i = 1, 2, . . . , 5 below (which will depend on ck, a, and 
k) is a Lyapunov covering for all a G and all k, G ]R>o x (0, ^2(0:)]. 

To allow further flexibility in Stage 1 and Stage 2 of the procedure, there 
is a multitude of parameters. We have thus provided a list of all parameters 
and their chosen values in the appendix. Moreover, to illustrate the regions 
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Ui ior i — 1,2,... 5, we have provided Figures 3.2-3.5. In what follows, the 
most crucial ingredient is where Ui covers the unstable trajectory 

along the negative X-axis. 

Lyapunov Function 1. Define 

M^,y) = Cii5\xf-u'\xf+'') 

on 

Ui^{x< -2} n {\y\ < 2\x\''^^^} . 
Then (pi is a Lyapunov function on Ui for all a e and all Ki > 0. 

Proof. The region excludes the y-axis, so 0i G C°°{Ui). Note by the choice 
of Ui, (j)i{x,y) > Cilxl'^, therefore (j)i{x,y) — )■ oo as {x,y) — )■ oo, {x,y) G Ui. 
Thus (I) and (II) are satisfied in t/i. To see (III), note that: 

L(t>i{x,y) < -Ci{2n2-haMxf^^ + 0{\xf). 

Therefore there exists > such that: 

L4>i{x,y) < -^(2«;2-5q;i^)|x|^+i + L'i (3.1) 
- -]^(2'^2-5ai/3)0i(x,y) + L'i 

where the last inequality follows by the choice (see the Appendix) of cr G 
(ai/a2, 1), 5 = ^, and /3 = (2 + a)5. □ 

Lyapunov Ftinction 2. Define 

<t>2{x,y) = C2 

on 

U2^{x< -2} n < \y\ < 2} . 

Then 4>2 is a Lyapunov function on U2 for all a e and all ki > 0. 



\x\ 


2<5 _|_ 


\y\ 


126 




\y\ 


\2aS 
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Proof. This region excludes both x and y-axes, therefore 02 ^ C°°{U2). 
Note that (l)2{x,y) > 4:~''^C2\x\'^^ on U2, whence ^2(3^, I/) 00 as {x,y) — >■ 
00, {x,y) G ?72. Thus 02 satisfies (I) and (II) on 1/2- To see (III), note that: 



L(t)2{x,y) = 2C25^|^((ai-aa2) + /t2fx(2a(5 + l)|x|-^|y|-2 + o(l)) 



1 2; 1 25+1 

< 2C25^|^ ((«i - cr«2) + /^2fT(2a(5 + 1) + 0(1)) (3.2) 
By the choice of K2, 

(ai — aa2) + K2cr{2a6 + 1) < 0. 
Therefore there exists D2 > such that 

L(f)2{x,y) < 6{{ai - aa2) + K2a{2a6 + l))(f)2{x,y) + D2. (3.3) 



Lyapunov Function 3. Define 

'Dx^ + y^Y 



03 (X,?/) = C3 



|2(T 



□ 
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on 



Us = {x<-l/2}n{\y\>l}. 
For all a e and all k & M."^ , (p^ is a Lyapunov function in U3. 

Proof. This region excludes the x and y-axes, so (63 G C°°{U3). Moreover, it 
is clear that 03(a;,y) — > 00 as {x,y) — )■ 00, {x,y) G U3. Thus 03 satisfies (I) 
and (II) in Us. Since 5 G (0, 1), we note that the S{S — 1) terms in dxx4>3{x, y) 
and dyy(j)3{x,y) are negative. Thus we have 



L(j)3{x,y) < 2d(j)s{x,y) 



{aix — aix"^ + y 



Dx 



+ (oay - a2xy) 



Dx^ + y^ 
aDsgn.{y)x^ + {I — a)sgn.{y)y'^ 



\y\{Dx^ + y^) 



K2a{2a + l)Dx^ . 
y^(Dx^ + y^) { 



< 2503(x,y) 



D{ai — aa2) 



\x\ 



Dx'^ + y'^ 



(3.4) 



+ {{D - a2{l - a))x + \a2\{l - a)) ^J^ 



+ 



0(1) 



By the choice of a G (q;i/q;2, 1) and D > max(l, {a2 + 2|a2|)(l — cr)) there 
exist constants C'^^D^ > such that 



L(t)3{x,y) < -C'^(t)3{x,y)+Ds, 
for {x,y) G t/3. This proves (III) is valid for ^3 in Us. 

Lyapunov Function 4. Define 

cj),{x,y) = C,{-x + N\y\'^'-'^^') 

on 

U, = {-l<x<N}n{\y\>l}. 
For all a, K & M^, 04 is a Lyapunov function in U4. 



□ 
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Figure 3.3: The regions U2 and U3 



U2 




U2 n U3 


- U2 — — -- 






_fi _4 _3 


2 -1 


U2 ■ 


U2 n f/3 


t/3 


U3 



Proof. This region excludes the x-axis, so 04 G C°^(f/4). As (a;,y) — ?■ 00 
with {x,y) G f/4, |?/| — t- 00 since a; is bounded. Therefore 04(x,?/) — )■ cxd as 
(x,?/) — )■ 00, G f/4. Thus 04 satisfies (I) and (II) in f/4. To see (III) is 

valid in f/4, note that 

Note this implies that there exist C4, Z^4 > such that 

L(j)4{x,y) < -C404(a;,y) + -D4, 

for {x,y) G f/4. 

□ 

Lyapunov Function 5. Define 

05 (a;, y) = C5(?7x^ + ?/^)'^ 

on 

f/5 = {a; > iV/2}. 
For all a, K G M^, 05 a Lyapunov function on f/5. 
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Figure 3.4: The regions U3 and U4, with > 1. 
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U3 



Us n U4 



-2.0 -1.5 



Us 



.0 -0.5 



Ui n U4 



U4 



0.5 1.0 



U4 



Proof. This region excludes the origin, so 05 G C'^{U^). Moreover 05 — )■ 00 
as (x,?/) — 7- 00. Thus 05 satisfies (I) and (II) in f/5. Since 7 G (0, 1) we see 
that the 7(7 — 1) terms in dxx(p5{x,y) and dyy(t)^{x,y) are negative. Thus we 
have 



L(j)5{x,y) < 2705(x,?/) 

+ {a2y - a2xy) 
< 2705(a;,?/) 



aix — a\x^ + y^ 



rjx 



y 



2 

aiTjx 



rjx^ + y^ 
3 



rfx"^ + y"^ 
+ 0(11 



r]x^ + 



((a2 - v)x - |a2|) 



?7x^ + y^ 



+ 



|ai|?7x 



+ 0(1) 



By the choice of = 0:2/2 and > 4|a2|/«2, there exist constants C5, D5 > 
such that 



L(j)5{x,y) < -Ci^(f)5{x,y) + D5, 
for (x,?/) G f/5. Hence 05(x,?/) satisfies (III) in t/5. 



□ 



By construction, the regions Ui, U2, U3, U4, f/5 overlap in such a way that 
we can extract a strong Lyapunov covering {(f)i,Vi), (02, V2), (05, V5) 
subordinate to the Lyapunov covering (0i,f/i), (0i,f/2), (05,^/5). We 
now proceed onto Stage 2. 
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Figure 3.5: The regions U4 and with N = 2. 
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3.3 Stage 2 

We now prove that 0i, 02? • • • ? 05 constructed above can be "patched" to- 
gether to yield a globally-defined Lyapunov function $ as in Lemma 2.3 By 
Proposition 3.2, to prove Lemma 2.3 it is enough to show that for all 02 > ai 
there exists k,2 = '^2(ck) > such that 

for all a G and all Ki > 0. In what follows, for all 0:2 > ai, we construct 
$ = such that for some ^2 = «^2(ck) > 0: 

*a,o e Ar(a,a,K), 

for all a G and for all ki > 0. Hence the choice of K2 can depend on a 
but not on a nor on Ki. 

To obtain $ as above, we define smooth auxiliary functions 

: Ui U Ui+, ^ [0, 1] 

such that Pi^i+i = on f/j and Pi^+i = 1 on f/j+i for i = 1, 2, . . . , 5, and 
prove that 



■= Pi,i+l4>i+l + (1 - Pi,i+l)4>i 



(3.5) 



17 



is a Lyapunov function on Ui U t/j+i for i = 1, 2, . . . , 5. Note then we can 
choose $ G C°°(]R^) such that outside of a sufficiently large ball Bn, we have 



on [/A (U,ti(f/, n « = 1, 2, . . . , 5 

-1 on Ui n i = 1,2, ... ,5. 



It follows then by construction, 

1. ^{z) — > oo as l^l — )■ oo. 

2. There exist constants C,D > such that L^{z) < -C^{z) + D for all 

z e M^. 



Note that by expression (3.5), for the newly defined we do not need 
to verify (I) and (II) in Definition 3.1 All we need to do is prove (III) on 
Ui n f/j+i. To simplify matters later, we note that for i = 1, 2, . . . , 5: 

L(f)i^i+i = pi^i+iL(f)i+i + (1 - pi,i+i)L(j)i (3.6) 
+ (aix - aix^ + y^){4>i+i - (pi)d^pi^i+i 
+ {a2y - oi2xy){(j)i+i - 4>i)dyPi^i+i 

+ Kl(0i+l - (j)i)dxxPi,i+l + 2Kidx{(j)i+l - (f)2)dxpi,i+l 

+ l^2{4>i+l - 4>i)dyyPi,i+l + 2K,2dy{(f)i + l " 02 ) <9j; Pj,i + 1 , 

on Ui n Ui+i. 

Patch 1. First we will patch together 0i and 02 to yield 0i^2, a Lyapunov 
function on Ui U U2- To do this, we first define the function pi^2- Let q : 
^ M and /, ^ : M ^ M be defined by 



q{x,y) = |x|^/^|y|-l, 

'exp(Trp^) if t e [0, 1] 
if [0,1] 

9it) = ^ I f{s)ds. 



m 



|1 J-oo 

Let pi,2 : U1UU2 ^ [0, 1] be defined by 

piA^^y) = 9{Q{x,y))- 

It is easy to see that pi 2 G C°°{Ui U U2) satisfies pi^2 = on f/i and pi^2 = 1 
on U2. 
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Claim 1. Let 0i 2 = Pi,202 + (1 ~ Pi,2)0i- There exists ^2(0:) > such that 
01,2 is a Lyapunov function on Ui U U2 for all a G and all Ki > 0. 



Proof. By estimates (3.1) and (3.3) and our clioice of parameters, we liave 

p^^2L(p2 + {l- Pi,2)L<pi < -/€2^i|x|(2+-)^+i + max(Di,D2), (3.7) 

for some constant Ei > independent of botli a and k. For i = 1 the 
remaining terms on the right-hand side of relation (3.6) are equal to 

-aiX^(02 - 0l)5xPl,2 - Oi2Xy{(f)2 - 0l)5j/Pl,2 (3.8) 
+ /t2(02 - h)dyyPl,2 + 2/^2^,(02 " h)dyPl,2 + Oi\x\^^+''^'+'/^) 

Since Ci > C2, we see that for sufficiently large {x,y) G f/i fl f/2: 
-aix2(02-0i)9.pi,2-«2a;i/(02-0i)9,pi,2 < -^2/(g(x, y))|x|(2+-)^+i 

(02 -0l)5>l2 + 29,(02 -0l)5,Pl2 < ^J^^l^l''^"''^' 

where h{t) = t^{t — 1)^ and E2,E3 > are independent of both a and k. 
Note that the quotient f{t)h(t)^^ — t- as t — ?■ or 1. Therefore there exists 
e > independent of a and k such that for Ixl""^/^!?/! G (1, 1 + e) U (2 — e, e) 
we have 

E,f{q{x,y)) ^ ^^^3^ 

This implies that for sufficiently large (x, y) E UiD U2 with la;!"*^/^!?/! G (1, 1 + 
e) U (2 - e, e), (Q is bounded by 

-K2Ei/2\x\^^+^^'+\ 

Note that l/h{q{x, y)) is bounded and f{q{x, y)) bounded away from zero on 
G [1 + e, 2 — e]. From this, we see that there exists ^2 > sufficiently 
small (which only depends on ex) for the estimate 

^2/2 > 



h{q{x,y)) 

to hold on G [1 + e, 2 — e]. These estimates imply that for sufficiently 

large (x, y) G t/i fl U2 with G [1 + e, 2 — e], (3.8) is bounded by 

-C"|x|(2+-)^+i 

for some C > 0. This finishes the proof of Claim 1. □ 
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For the remainder of Stage 2, let A; : M — )■ M be a smooth function such 
that k{t) = for t < 0, k{t) = 1 for t > 1, and k'{t) > on (0, 1). 

Patch 2. We will now patch 02 and 03 together to obtain a Lyapunov func- 
tion on U2UU3. To this end, let p2,3{x, y) = k{\y\ — 1). It is clear that k can 
be chosen so that p2,3 e C°°(f/i U U2). 

Claim 2. Lei 02,3 = P2,303 + (1 ~P2,3)02- Then for all a G and all ki > 0, 

02.3 zs a Lyapunov function on f/2 U f/3. 

Proof. Since p2,3 is independent of x, 

L4>2,3 = P2,3L4>3 + {1 - P2,3)L4>2 + {0'2y - a2Xy){(j)3- (f)2)dyP2,3 
+ /t2(03 - 4>2)dyyP2,3 + 2K2dy{(t>3 - 02)<9j/P2,3- 

We chose C3 < €2/0^ so that, for sufficiently large |a;| with (x, y) G t/2 H t/3, 

-a2X|/(03 - (t)2)dyp2,3 < 0. 

Note moreover that, by (3.2) and ( |3.4 ), on Ui fl f/2, we have 

^02,3 < - a2a;y(03 - 02)9,P2,3 + 0{\x\^'), 

which implies the claim. □ 

Patch 3. We now patch together 03 and 04 to obtain a Lyapunov function 
on f/3 U f/4. Define p3,4 : ^ M by P3,4(x, y) = k{2x + 2). 

Claims, i^ei 03,4 = p3,404 + (1 — p3,4)03. Then for all a E M.'^ and all Ki > 0, 

03.4 is a Lyapunov function on f/3 U f/4. 

Proof. Since p3,4 is independent of y, 

L(j)3A = P3AL(f)4^ + - P3A)L(j)3 + {aix - aix^ + y^){(f)4- (p3)dxP3,4 

+ ft:i(04 - (p3)dxxP3A + 2fi:i9^.(04 - (p3)dxp3A- 

By the choice of NC4 < C3, 04 — 03 is negative for sufficiently large y in the 
region. Therefore, there exist E^,Eq,Ej > such that 

L(t)3A < -p3,4i55y'-((l-P3,4)i?6-i57|5x.x.p3,4|y'^'-'^^ + C^(l).(3.9) 
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Since fc, fc" — )■ as t — )■ there exists e > such that 

(1 - P3,4)^6 - ^7|<9xxp3,4| > ^8 

for all {x,y) G U3 H with x G (—1,-1 + e) and some > 0. For 
{x,y) G t/3 n f/4 with x G [—1 + e, —1/2), p3_4 > i^g for some Eg > and 



|c^xxP3,4| is bounded. Putting these two estimates together with (3.9) finishes 



the proof of the claim. □ 
Patch 4. Lastly, we patch together (j)^ and 05 to obtain a Lyapunov function 



64^5 on f/4 U f/5. Define then p4^5 : — )• M by Pi^^{x,y) = k{jfX — 1). 



Claim 4. i^ei 04^5 = p4_5</)5 + (1 — p4_5)(/)4. Then for all a E M."^ and all ni > 0, 
04,5 zs a Lyapunov function on f/4 U f/5. 

Proof. With the choice of NC^ > C5 the proof is nearly identical to the proof 
of Claim 3, so we omit further discussion. □ 

4 Regularity and Controllability 

We now prove Lemma [2^ We first start by showing that the transition mea- 
sures P{t, z, ■ ) have smooth densities with respect to Lebesgue measure on 
M^. This is a simple consequence of Hormander's hypoellipticity theorem [o], 
a powerful result that says, roughly speaking, that if enough noise spreads 
throughout the system such regularity follows. As we shall see, the word 
"enough" can be quantified by considering the span of the Lie algebra gen- 
erated by vector fields in the expression for the generator L. In general, this 



requires rewriting the SDE in the Stratonovich form 18 . In the present case 



Ito and Stratonovich forms coincide, since the noise coefficients are constant. 



4.1 Regularity 

The regularity of the measures P{t,z, ■ ) can be studied using either classical 



PDE theory [9] or the Malliavin calculus 18 19 . To state the result, for 
vector fields X and Y, let f-'^, ^] be their Lie bracket (commutator). If we 
write the generator L in the form: 

L = Xq + Xi + X2, 
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where Xq, Xi, and X2 are vector fields (treated as first-order differential 
operators, so that X| denotes a composition of such operator with itself) 
with smooth coefficients and if at each point 2; G the list: 

Xniz) Ji = l,2 

[%,^,.](^) Ji,j2 = 0,1,2 

[X^„[X,„X,^]]{z) ji,j2,j3 = 0,1,2 



spans M^, then for t > 0, ^ G 

P(t, z, dw) = pit, z, w) dw, 

where p{t, z, w) is a smooth function for (t, z, w) G (0, 00) x x M^. More- 
over, if u is an invariant probability measure, then 

p^dw) = Pviw) dw 

where G C°°(M2). 

Proof of Lemma 2.4\ (!)■ Note that 



d d 



Thus let 



f d d 

Xq = \{aix - aix'^ + y^)— + {a2y - a2xy) — 

X2 = 



'dy' 

If both Ki > 0, K2 > 0, then Xi{z) and X2{z) span for all z G M?. If 
Ki = 0, we see that 

d 

^2,2,0 '■= [-^2, [X2,Xo]] = 2^2 — . 
Since K2 > 0, X2,2,o(-2) and X2{z) span for all z G M^. □ 
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4.2 Controllability 



We now prove Lemma 2.4 (2) by utilizing a connection between SDEs and 
control theory provided by the Stroock-Varadhan support theorem 22 . We 
will also use and modify results in [l , 11, 12 to work in our setting. 

With fixed ai,a2 ^ 0^2 > ai > 0, and ki > 0, ^2 > 0, consider the 
ordinary differential equation: 



xt = {aixt - aixl + yl) + Kiu{t) 

yt = {a'2yt- oi2Xtyt) + i^2v{t), 



(4.1) 



where m, v : [0, 00) — t- M are called controls. We assume u and v belong 
to W, the class of piecewise constant mappings from [0, 00) into M with at 
most finitely many discontinuities. For w{t) := (^(t), v{t)) with u, v ^U, we 



let z{zo,'w,t) be the maximal right integral curve of (4.1) passing through 



at t = 0. 

In view of |l], one can use accessibility sets of the system (4.1) to help 
determine the support of an extremal invariant probability measure. More 
precisely, for t > 0, we let A{zq, < t) denote the accessibility set from zq in t 
units of time or less, i.e., 

< t) := {z{zo, w,s) : w = {u, v) , u, v G U , < s < t}. 

Defining 

A{zo) :=\jA{zo,<t), 

t>o 

we note the following result in [ij: 

Theorem 4.1. If fJ^ is an extremal invariant probability measure for the pro- 
cess Zt, there exists zq G such that 



supp fi = ^(2:0). 

Proof. This follows by Proposition 1.1 and Definition 1.1 in |[l] since Zt has 
continuous density with respect to Lebesgue measure. □ 



Using Proposition 2.1 of [I] and Lemma 2.4 (1), to prove uniqueness of 
the invariant probability measure ly it is enough to show: 
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Lemma 4.2. For all Zq,Wo G M^, there exists a non-empty open subset U = 
U{zo,Wo) C such that: 



A{zo)nA{wo) D U{zo,wo) 



(4.2) 



To prove the lemma above, we require ideas in (TT, 12 . This is because 
it is difficult to find specific controls u, v and then solve equation (4.1) to 



determine A{zq). As emphasized in these works, one can simplify such a 



procedure for systems like (4.1) by using geometric control theory. 



To illustrate this approach, let Z and W be the vector fields on 
mined by 



deter- 



Z{x,y) 
Wi{x,y) 
W2{x,y) 



{aix — OL\x^ + ?/^, a2y 



a2xy) 



(/€l,0) 
(0,^2), 



and define a family of smooth vector fields: 

F = {Z + uWi + vW2 : u, i; G M} . 



2 let exp{tZ){zo) 
= 0. For zo G 



We call F a polysystem. If Z G -F, for t > and zq E ^ 
denote its maximal integral curve passing through zq at t 
and t > we let Ap{zQ, < t) denote the set of G such that there exist 

Zi, ^2, . . . , Zk & F and times ti,t2, . . . ,tk > such that ti + ^2 H \-tk <:t 

and 

expitkZk) o exp(tk-iZk-i) o ■ ■ ■ o exp{tiZi){zo) = w. 

It is easy to see that Ap{zo, <t) = A{zq, < t) for all zq G M^, t > 0. 

The benefit of using geometric ideas is that it provides a means by which 
to enlarge F without changing Ap{zQ, < t). In light of this, we say two 
polysystems Fi and F2 are equivalent, denoted by Fi ~ F2, if for all zo G 
and t > 

AFAzo,<t) = Ap^{zo,<t). 



One can show, see 11, 12 , that if Fi ~ F and F2 ~ F, then Fi U F2 ~ F. 



Thus to make F as large as possible, we consider the union of all polysystems 
equivalent to F, called the saturate of F, which we denote by Sat(F). To yield 
new equivalent polysystems from old, we require the following definitions and 



lemmata 11 12 . 
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Definition 4.1. We call a diffeomorphism : — t- a normalizer of the 
polysystem F if for all zq G and all t > 0: 



r]{AF{v-'{zo),<t)) C AF{zo,<t). 
We denote the set of all normalizers of F by Norm(F). 
Lemma 4.3. If F is polysystem, then 

U {v.{V) : V e F} r-. F, 

TjaNorm^F) 

where rj^, denotes the differential. 

Lemma 4.4. If F is a smooth polysystem, then F is equivalent to the closed 
convex hull of {XV : < A < 1, \^ G F}. Here the closure is taken in the 
topology of uniform convergence with all derivatives on compact subsets of 



With these ideas in place, we now prove Lemma 2.4 (2) 



Proof of Lemma \K4\ (2). Set F = {Z + UW1+VW2 : u,v e U}. By Lemma 



4.4, for all constants A G 



XWi = lim -{Z + nXWi) G Sat(F) 

n— >oo n 

XW2 = lim -{Z + nXW2) G Sat(F). 

It is easy to see that z 1— )■ exp{XWi){z) G Norm(F) for A G M, i = 1, 2. Since 
K2 > Lemma 4.3| implies: 



exp(Aiy2)*(^) = Z + X[W2,Z] + ^ [W2 [W2, Z]] G Sat(F). 



for all A G M. Applying Lemma 4.4 with /i G M, we see that 



lim hexp{X\2fi\'/^W2UZ)) = [W2 [W2, Z]] G Sat(F). 

A— >-oo A 

Moreover, [W2 [W2, Z]] (2) = (|/i|K2,0). Using this vector field with W2, 
we see that for zq G and t > 0: 



AF{zo,<t) D H{zo), 
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where for z = {x,y), H{z 
proof of Lemma 



4.2 



{w = (-Ujf) G : u > x}. This finishes the 



(2) since for all zq, ujq G 



A{zo) n A{wo) D U {zo, Wo) ^ 0, 
for some U{zq,wq) open. 



□ 



5 Instabilities 



We assume now that ai > a2 and prove Theorem 2^ Similar to Section |3} 



this will be established by constructing an appropriate test function \E', as in 



the hypotheses of Lemma 2.5 Before we proceed on to the construction of \E', 



we first discuss the deterministic dynamics under the assumption ai > a2- 
5.1 A Robust Explosive Region 

It is not hard to see that when a2 > Oii as in Section [3} the deterministic 
dynamics 

x{t) = aixit) - aixitf + yitf (5.1) 
y{t) = a2y{t) - a2x{t)y{t) 

has a single unstable trajectory along the negative x-axis. When ai > a2, we 



prove this instability in equation (5.1) is more robust. Remaining consistent 



with the earlier notation, we let z{t) = {x{t),y{t)) be the solution to equation 



(5.1 ), defined on the maximal interval of times t > on which it exists. The 
initial condition will always be made clear in the following arguments. We 
first need a definition. 

Definition 5.1. We call a set f/ C invariant if for all z{0) = {xq, yo) G U, 
z{t) eU for all times t > for which the solution z{t) is defined. 

Using Lyapunov-type criteria, we will determine an invariant set C 
for which solutions starting at zq G A^ escape to infinity in finite time. To 



this end, we recall that the infinitesimal operator for (5.1) has the form: 



d d 
LP = (aix - aix'^ '^y^'^dx^ ^^^^ ~ ^'^^^'^dy' 
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For ^, M > 0, define functions V'l, V's : by: 

i^2ix,y) = iy-x 
i^3{x,y) = -^y-x, 

and let 

u^,M = i'l'm, 00)) n V2~'((o, 00)) n ^^\io, oo)). 

Proposition 5.1. For all ^ > such that < cti — a2, there exists M > 
such that U^,M is invariant. 

Proof. Since < < ai — 0:2, we may choose M > sucli tliat for all 
i = 1,2,3: 

L%{x,y)>0, 
for all {x,y) e U^,m- Fix 2;(0) = (a;o,2/o) G f^5,M and let 

S = inf{z(t) e \ U^,m}. 

We have for ti < ^2, ^i, ^2 e [0, -S"): 

-X{t2) + X{ti) = ^Iji{z{t2)) - = / i^Vl('2H)c^M > 0. 

Thus x{t) is strictly decreasing on [0,5'). By continuity, z{t) cannot exit 
U^^M through the vertical line x = M. Using similar reasoning, we see that 
4'2{z(t)) and ips^zit)) are strictly increasing on [0,5'). This implies that z{t) 
cannot exit U^^m through the lines \x\ = ^\y\. □ 

For ^ > such that < ai — a2, choose M > so that the conclusion 
of the preceding proposition is satisfied and 

LVi(x,y)>C[V'i(x,y)]2, (5.2) 

for all {x,y) e U^,m for some C > 0. Let ^4° — U^^m- We now have the 
following proposition, illustrated by Figure 5.6: 

Proposition 5.2. For all z{0) = (a;o,|/o) G A^, z{t) reaches infinity in finite 
time T > where 

T< 1 



C\xq\ 
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Figure 5.6: The region ?7i^i5 and approximate dynamics when ai > a2- 
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Proof. Fix z{0) = {xQ,yo) e A° and let T„ = mit^Q{z{t) G -B^} where Bn is 
the open ball of radius n centered at the origin. Let T be the finite or infinite 
limit of T„ as n — )■ oo. By the previous proposition, z(t) E = U^,m for all 
t E [0,T). Fix t G [0,T). By the proof of the previous proposition, the map 
s H- x{s) is strictly decreasing and smooth. As such, it is a diffeomorphism 
onto its image. Using to denote its smooth inverse, we have: 



< 



x-\x{t)) - x-\x{0)) 

-^x~\y) dy 
x(0) "Z/ 

x{t) ^ 
x{0) ^ 
x{0) ^ 

dy 



it) C{x{x~\y)y- 



1 

c 



\xo\ \x{t)\^ 



(5.3) 



< 



C|xo| 
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where (5.3) follows from (5.2). Since t G [0,T) was arbitrary, we see that 

1 



T < , , 
C\xo\ 



□ 



5.2 The Random Dynamics 

The robustness of the explosive region U^^m suggests that noise will not 
be sufficient to stabilize the system, a fact which we will now prove. Let 
q : (-00, 0) X M ^ M, / : M ^ M, p : (-oo, 0) x M ^ M, and /i : ^ M by 

/ ^Mt^) ifte [-1,1] 



<l{x,y) = —, f{t) 



X 



ift^[-l, 1] 



r 1 7 / \ I exp(— It^) if X < —M 



if X > -M ■ 

Finally, let g = hp (the product of h and p) and note that g is smooth on 
M."^, non-negative, bounded, and is supported in U^^m- Moreover, g is strictly 
positive on U^^m- We will use g to show that the system in not regularized 
by noise in the specified parameter range. 

Proposition 5.3. For all > such that {ai — a2)/2 < < — «2 there 
exist constants M, C > such that Lg > Cg on M? and g is strictly positive 
on U^^M- If o-i > 0-2 then M can he chosen independent of ^. 

Proof. First note that the function h only depends on x; hence, applying L 
to g we obtain: 

Lg = pLh + {aix — aix'^ + y'^)hdxP + {a2y — a2xy)hdyp + KihdxxP 
+ K,2hdyyp + 2Kidxpdxh. 

To estimate Lg, we find it convenient to use the following expressions for the 
derivatives above: 

aix'^ — aix — y"^ 2(x + M) + l,, 
' (x + M)2 ^ (x + M)4 ' ' 



-1 2g2 -2g2 

2g2(g4 + 4^2 _3) 2e{3q*-l) 
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all of which follow from the definitions of h and q. Let Di > and D2 > 
be upper bounds for 

2t2(t4 + 4^2 _ 3) 2(3i^ - 1) 

and 



(1-^2)4 (l-f2)4 

respectively for t G (—1, 1). Therefore, since dxpdxh > and > (^q;-^ _ 
q;2)/2, we have the following estimate for the second-order terms 

i^ihdxxP + K2hdyyp + 2Kidxpdxh > -{kiDi + K2D2^'^)x~'^g 

> - KiDi H X g, 

\ ax-a2j 

on U^^M- To handle the remaining terms we first recall that y'^x~'^ < /^"^ on 
t/^^M- Hence, 

{aix - aix'^ + y^)hdxP + (02^ - a2xy)hdyp > 9 ^-^ _ ^2^2 ( ~ "2 

-C"^ + kr^(ai - as)^, 

and 

r, ^ /(ai-r^)a;2-aix 2Avi(a; + M) + 1\ 

^^^^l — i^TJIy — + ix + Mr ) 

Let D3 > he an upper bound for 

2ki{x + M) + 1 
(x + M)2 

on (— 00,— M). Therefore 



^2 



^ [x + My 



Together these estimates imply: 



(a2 - \a,\M-' - D^M-^ - ( t^rD, + ^^^^'j m-\1 + M/xf , . ^ . 
\ \ ai- a2j J {x + My 

2\x\q'^(ai - Q-2- + |x|"^(ai - 02))" 



[l-q 



2\2 
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For all > 1 and M sufficiently large (the choice of which does not depend 
on 0: 

T y f «2 2|x|g2(ai - ^2 - ^"^ + (qi - a2)\x\~'^) \ 

^ -\r]{l + M/xy (l-g2)2 J 9 

for all C, satisfying ai— 0:2 > > (ai — a;2)/2. If ai > a2 the second term on 
the right is nonnegative; thus we have proved the proposition when ai > 02- 
If 02 > cti, we have 

for some constant C^/ > 0. For a given satisfying ^"^ < — ^2 the second 
term is positive for sufficiently large M, thereby proving the proposition in 
this case as well. □ 



6 Conclusions 



We saw that in equation (1.1), ai = 02 represents a critical barrier for 
ergodicity. In particular, it was shown that if 02 > ai, there exists a unique 
invariant probability measure. This result was essentially proven by the 
existence of a globally-defined Lyapunov function $ such that L$ satisfies 



the strong bound (2.2). The function $ was constructed by first exhibiting 



a strong Lyapunov covering {{(pi, the existence of which, at the very 

least, guarantees that Zt cannot exit directly through any Vj. To eliminate 
the (highly unlikely) possibility of other explosion scenarios, the covering was 
patched together to yield $. 

In the case when 02 > «i and ki > 0, we were able to quantify the rate of 
convergence to the steady state. As we saw, convergence is exponentially fast 
in a norm stronger than the total variation norm on i3-measures. If Ki = 0, 
ergodicity still remains true in this case, but it is much harder to quantify 
the rate of convergence. This should be expected, since there is less noise in 
the system. Methods developed in [T||9|[Tl|[T2|[22] proved useful in showing 
ergodicity in this situation. 

In the case when ai > a2, we saw that ergodicity cannot even be dis- 
cussed, since global stability is not satisfied. This is because solutions start- 
ing in a wedge-like region containing the negative X-axis explode in finite 
time with positive probability. We believe that this is primarily due to the 
presence of a more robust explosive region in the deterministic dynamics. 
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8 Appendix 



In this section, we list the parameters introduced in Stage 1 of Section [3] 
along with their chosen value. Recall that Q!2 > cti > are fixed constants. 

a = — G(ai/a2,l) 

2^2 



6 



_K2_ 

8ai 



fj = {2 + a)d = —- K2 

Lbaia2 

, NT (02 -ai) 

7 = (l-^)^ = ^7^ 

1dci;iQ;2 



a2 

y 

(02 + 2|a2|)(a2 - ai) ^ . , ^1 .w-, ^^ 

D = IH > max(l, (0:2 + 2|a2|)(l — o-j) 

a2 

4|a2l , , , , X 

= 1 + — ^ > max(l,4|a2|/a2) 
a2 

Ci = 2 
C2 = 1 

C4 
C5 



2D^ 
1 

1 



With these choices, choose ^2 = ^2(0!) > such that: 

5, 7 e (0,1/2), 

(ai — (Ta2) + K2a(2(y5 + 1) < 0, 

< E2/2 for all \x\^^^\y\ G [1 + e, 2 - e] 



h{q{x,y)) 

where e, E2, E3 > and h{q{x, y)) were introduced in Patch 1 of Stage 2 and 
are all independent of k and a. It is easy to see that K2 only depends on a 



which is permissible by Proposition 3.2 
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